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Abstract 

We consider a noncommutative version of the Davey-Stewartson equations and derive two 
families of quasideterminant solution via Darboux and binary Darboux transformations. These 
solutions can be verified by direct substitution. We then calculate the dromion solutions of the 
equations and obtain computer plots in a noncommutative setting. 

1 Introduction 

The Davey-Stewartson (DS) equations have become a topic of much interest in recent years. De- 
rived by A. Davey and K. Stewartson in 1974 [3J, the system is nonlinear in (2 + l)-dimensions 
and describes the evolution of a three-dimensional wave-packet on water of finite depth. By car- 
rying out a suitable dimensional reduction, the system can be reduced to the (1 + l)-dimensional 
nonlinear Schrodinger (NLS) equation. 

A major development in the understanding of the DS equations came in 1988, when Boiti et 
al. [3] discovered a class of exponentially localised solutions (two-dimensional solitons) which 
undergo a phase shift and possible amplitude change on interaction with other solitons. These 
were later termed dromions by Fokas and Santini [6], derived from the Greek dromos meaning 
tracks, to highlight that the dromions lie at the intersection of perpendicular track-like plane 
waves. 

Multidromion solutions to the DS system have been obtained using a variety of approaches - the 
inverse scattering method [6], Hirota's direct method [13] and others. These solutions have been 
determined both in terms of Wronskian [13] and Grammian [10] determinants. 
Additionally, there has been considerable interest in noncommutative versions of the DS equa- 
tions. Hamanaka [12] derived a system with noncommutativity defined in terms of the Moyal 
star product [19] . while more recently, Dimakis and Miiller-Hoissen [5] determined a similar sys- 
tem from a multicomponent KP hierarchy. This then enabled calculation of dromion solutions 
in the matrix case. 
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The strategy that we employ here, whereby we introduce noncommutativity into an integrable 
nonlinear wave equation without destroying the solvability, has previously been considered by 
others in the field, for example by Lechtenfeld and Popov [16], and by Lechtenfeld, Popov et al 
in [15], where a noncommutative version of the sine-Gordon equation is discussed. 
In this paper we are not concerned with the nature of the noncommutativity, and derive a sys- 
tem of noncommutative DS equations in the most general way by utilising the same Lax pair 
as in the commutative case but assuming no commutativity of the dependent variables. This 
method has also been employed by Gilson and Nimmo in [11] for the case of the noncommutative 
Kadomtsev-Petviashvili (KP) equation. We find that the noncommutative DS system obtained 
in this manner corresponds to that given in an earlier paper by Schultz, Ablowitz and Bar Yaa- 
cov |21j . where a quantum version of the DS equation is ultimately discussed. 
We derive quasiwronskian and quasigrammian solutions of this system in section [4] via Darboux 
and binary Darboux transformations and, in section [U verify these solutions by direct substitu- 
tion. 

We then use the quasigrammian solution to determine a class of dromion solution and, by 
specifying that certain parameters in the solution are of matrix rather than scalar form, ob- 
tain dromion solutions in the noncommutative case. We conclude with computer plots of these 
dromion solutions. 



2 Noncommutative Davey-Stewartson equations 

We consider the system of commutative DS equations given by Ablowitz and Schultz in [2], with 
Lax pair 

L = d x - A + aJdy, (2.1a) 

M = d t -A + -Ad y - Udyy, (2.1b) 

where 

J-( l °) A-( ° ^y^) (2 2) 

J \0 -l)' A \r(x,y,t) J {Z - Z > 

for r = ±q* (q* denotes the complex conjugate of q) and A is a 2 x 2 matrix given by 



A 



M ^2 fa* ~ a %) 



i 



V 2cj 2 



(r x + ar y ) A 



(2.3) 
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We choose a = — 1 or a = i for the DSI and DSII equations respectively. By considering the 
same Lax pair (|2.1al b) as is used in the commutative case and assuming no commutativity of 
variables (we do not specify the nature of the noncommutativity), we obtain the compatibility 
condition 

-A x + [A, A] + A t - a J A y + -AA y - UA yy = 0, (2.4a) 
-A x + a[A, J] - UA y = 0, (2.4b) 



2 



from which we generate a system of noncommutative Davey-Stewartson (ncDS) eqations 
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m = ( qxx + a q w) + i( Aiq ~ ^2), (2.5a) 
in = ^ ( rxx + 0-2 r vv) ~ [ ( rAl ~ ^ 2r )' (2.5b) 



(d x + ad y )Ai = -— (d x - ad y )(qr), (2.5c) 
(8 X - ad y )A 2 = ^(8 X + ad y ){rq). (2.5d) 
(Note that we obtain from the above system the nonlinear Schrodinger (NLS) equation pQ 

i<?t + Qyy ± 2qrq = (2.6) 

and its corresponding complex conjugate by taking a dimensional reduction d x = 0, with 
A x = ±±qr, A 2 = ±\rq [H]). 



For notational convenience, and to avoid the use of identities later when verifying solutions, 
we introduce a 2 x 2 matrix S = (sij) = 1,2) such that A = [J, aS] [T7], and hence 




S= | 

Additionally, by setting 

A = —S x - iJS y , (2.8) 
equation (|2.4bp is automatically satisfied, and (|2.4ap becomes 

S X x + iJSS x — iSJS x — iS x JS + iS x SJ + \aJS v JS 

a y 

— iaJSySJ + aJSt — crStJ — iaJSSyJ + iaSJS y J + iaJS yy J = 0. (2.9) 

Note that this is essentially the noncommutative analogue of the Hirota bilinear form (see for 
example [14 ) of the DS equations. 



3 Quasideterminants 



Here we briefly recall some of the properties of quasideterminants. A more detailed analysis can 
be found in the original papers 0[9]. 

The notion of a quasideterminant was first introduced by Gelfand and Retakh in [9] as a straight- 
forward way to define the determinant of a matrix with noncommutative entries. Many equiva- 
lent definitions of quasideterminants exist, one such being a recursive definition involving inverse 
minors. Let A = (aij) be an n x n matrix with entries over a usually noncommutative ring 1Z. 

nnooirlofnT'mmovTf Vwr I /II.. 



We denote the (i,j) th quasideterminant by \A\ij, where 



\A\ VJ = a ij -rl{A^)- 1 s i j . (3.1) 
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Here, A lJ is the (n — 1) X (n — 1) minor matrix obtained from A by deleting the z th row and the 
j th column (note that this matrix must be invertible), r\ is the row vector obtained from the i th 
row of A by deleting the j th entry, and s*- is the column vector obtained from the j th column of 
A by deleting the i th entry. 

A common notation employed when discussing quasideterminants is to 'box' the expansion ele- 
ment, i.e. we write 

an 



\A 



11- 



(3.2) 

«21 CL22 

to denote the (1, l) th quasideterminant. It should be noted that the above expansion formula is 
also valid in the case of block matrices, provided the matrix to be inverted is square; for example 
considering a block matrix 

'N 5 N 
C d / 

where N is a square matrix over 1Z, B and C are column and row vectors over 1Z of compatible 
lengths, and d G 7Z, we have 

N B 

(3.3) 



C 



d 



d-CN^B. 



Quasideterminants also provide a useful formula for the inverse of a matrix: for an invertible 
n x n matrix A = (aij) (i, j = 1, . . . , n), the (i, j) th entry of A~ l is given by 



(3.4) 



When the elements of A commute, the quasideterminant | A is not simply the determinant of 
A, but rather a ratio of determinants: it is well-known that, for A invertible, the (j, i) th entry 
of A^ 1 is 



,-det A 13 



detA 



Then, by (|3.4p . we can easily see that 



\A 



,• detA 



det AV 



(3.5) 



in the commutative case. 



4 Quasideterminant solutions via Darboux transfor- 
mations 

4.1 Darboux transformations 

Here we give a brief overview of Darboux transformations. Further information can be found in, 
for example, [18] . 

We follow the notation given in [11] . Let 9\ , . . . , 6 n be a particular set of eigenfunctions of an 
operator L, and define = (0\, . . . , 6 n ) and = (of 1 M for i, j = 1, . . . ,n, the n x n Wronskian 
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matrix of 9±, . . . , 9 n , where ^ denotes the k th y-derivative. 

To iterate the Darboux transformation, let 9^ = 9± and <j)M = <p be a general eigenfunction of 
L[l] = L, with L[ X ] covariant under the action of the Darboux transformation Gg^ = d y — 9^9Z^. 
Then the general eigenfunctions 4>[2] f° r ^[2] = Gq^L^Gq^ are given by 



with 



(<%]) 
#[2] = 0[2] U— 02- 



(4.1) 
(4.2) 



Continuing this process, after n iterations (n > 1), the ra th Darboux transformation of is given 

by 



where 



(4.3) 
(4.4) 



4.2 Quasiwronskian solution of ncDS using Darboux transfor- 
mations 

We now determine the effect of the Darboux transformation Gg = d y — 9 y 9~ l on the Lax operator 
L given by (|2. la|) . with 9 an eigenfunction of L. Corresponding results hold for the operator M 
given by (|2.1b[) . L is covariant with respect to the Darboux transformation, and, by supposing 
that L is transformed to a new operator L, say, we calculate that the effect of the Darboux 
transformation L = GgLGg is such that 

A = A-a[J,9 y 9~ 1 ]. (4.5) 

Recalling that A = [J,aS], we have S = S — 9 y 9 _1 , and hence, after n repeated Darboux 
transformations, 

S[ n +l] = S [n]~ (0[n])y0[n] ' ( 4 - 6 ) 

where Sm = S, 0m = 9. We express SL+i] m quasideterminant form as 



S 



n+l] 



S + 













9! .. 


9 n 








(n-2) 


fl (n-2) 










• t?n 








(n-1) 


fl (n-l) 


1 












1 


9^ . . 


n(n) 










v n 









(4.7) 



where ^ denotes the fc th y-derivative. It should be noted here that each 9i (i = 1, . . . , n) is not 
a single entry but a 2 x 2 matrix (since the 9{ are eigenfunctions of L,M). The Wronskian-like 
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quasideterminant in (|4.7p is termed a quasiwronskian, see , and [7] for details of Wronskian 
determinants. 

For ease of notation, for integers i, j = 1, . . . , n, we denote by Q(i,j) the quasideterminant [11] 





0(i+t) 



/; 


e 3 















(4.8) 



where, as before, G = (0^ 1 )i ; j=i,...,n is the nxn Wronskian matrix of #i, . . . , and ^ denotes 
the fc th y-derivative, Q is the row vector (6\, . . . , 9 n ) of length n, and fj and are 2n x 1 column 
vectors with a 1 in the (2n — 2j — l) th and (2n — 2j) th row respectively and zeros elsewhere. 
Again each B{ is a 2 x 2 matrix. In this definition of Q(i,j), we allow i,j to take any integer 
values subject to the convention that if either 2n — 2j or 2n — 2j — 1 lies outside the range 
1,2, ... , 2n, then = fj = and so Q(i,j) = 0. Hence (|4.7|) can be written as 



S = S + Q(0,0), 



(4.9) 



where So is any given solution of the ncDS equations. Here we choose the vacuum solution 
Sq = for simplicity. 

It will be useful to express the quasiwronskian solution (|4.9p in terms of the variables q and r, 
the variables in which the ncDS equations ()2.5at -d) are expressed. We have 



5 = Q(0,0), 



(4.10) 



which gives, by applying the quasideterminant expansion formula (|3.ip and expressing each 9i 
(i = 1, . . . , n) as an appropriate 2x2 matrix 



92i-l P2i 



(4.11) 



for 



(x,y,t), ij) = ip(x,y,t), an expression for S in terms of quasiwronskians, namely 



S 



/ 


G 


k 




G 


eo 


\ 











cj){n) 









Q 


/o 




G 


eo 




V 









^(») 





/ 



(4.12) 



where 4>^ n \ denote the row vectors \ ■ ■ ■ , </ > 2n' > ) > (^l"^' • • • ' respectively. By com- 

paring with (12. 7p . we immediately see that g, r can be expressed as quasiwronskians, namely 



2a 



G e 



-2(7 



e fo 



(4.13) 
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4.3 Binary Darboux transformations 

In order to define a binary Darboux transformation, we consider the adjoint Lax pair of the 
ncDS system (I2.5al -d). The notion of adjoint can be easily extended from the well-known matrix 
situation to any ring 1Z. An element a £ 1Z has adjoint a', where the adjoint has the following 
properties: if d is a derivative acting on 1Z, = —d, and for any product AB of elements of, or 
operators on 1Z, (AB)t = St At. Thus the adjoint Lax pair for the ncDS system is given by 

Lt = -d x - vlt - - Jd y , (4.14a) 
o 

Aft = -Q t - At + kr(4 + At^) + iJd ra . (4.14b) 

We construct a binary Darboux transformation in the usual manner (see for example [18]) by 
introducing a potential tp) satisfying the relations 

n(0,^) tf = V f 0» (4.15a) 
n(<f>,i/}) x = -oifi J<f>, (4.15b) 
0(0,^ = 1(^^-^0), (4.15c) 

with ^ an eigenfunction of L, M and ^ an eigenfunction of L<, M>. This definition is also valid 
for non-scalar eigenfunctions: if $ is an n-vector and ^ an m-vector, then 0( < I ) , 1 I / ) is an m x n 
matrix. We then define a binary Darboux transformation Gq <p by 

G e , P = l-0Sl(9, P r l dy l P ^ (4.16) 

for eigenfunctions 6 of L, M and p of Lt, Aft, so that 

0[2] = G^,^, (</>[!]) = 0[i] -%]O(%],p[i]) -1 O(p[i],0[i]) (4.17) 

and 

^[2] = Ge [1]lP[1] (^[i]) = V[i] -P[i]^(%],P[i])~ t ^(%] ) '0[i]) t ) (4-18) 

with 

8[2] = 0[2]U-e 2 , P[2] = ^[2] [V-P2. (4-19) 

After n > 1 iterations, the ra th binary Darboux transformation is given by 

4>[n+l] = 4>[n] ~ 0[n]O(0[ n ],P[n]) _1 O(p[ n ], <£[„]) (4.20) 

and 

1p[ n+1 ] = V[n] - Plnft(Q[n], P[n})~^(0 [n] , ^[ n ])\ (4.21) 

with 

#[n] = 4>[n]\cf>^6 n , P[ n \ = ^[n]\^p„. (4.22) 
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4.4 Quasigrammian solution of ncDS using binary Darboux trans- 
formations 

We now determine the effect of the binary Darboux transformation Ge tP on the operator L given 
by (|2.1a[) . with 9 an eigenfunction of L and p an eigenfunction of . Corresponding results 
hold for the operator M given by (|2. lbj) and its corresponding adjoint M<. The operator L is 
transformed to a new operator L, say, where 



L = Ge, p LG 



-i 



(4.23) 



We find that A = A - a[J, 9 y 9~ 1 ] - <j[9y9~ 1 , J], and hence, since A = [J,aS], it follows that 

S = S-6Sl(6,p)- 1 pi (4.24) 

with 9 = —9Q(9, p)" 1 . After n repeated applications of the binary Darboux transformation Gq jP , 
we obtain 



S[n+1] - S[n] ~ ^2 0[k]M(9[k]i P[k]) V[fc]> 



(4.25) 



fe=i 



where Sm = S, 9m = 9. Defining 6 
quasigrammian form [11] as 



h,-..,9 n ) and P = (px, . . . , p n ), we express S[ n+1 ] in 

o(e,p) pt 



Si 



n+l] 



s 



e 
















(4.26) 



where is the Grammian-like matrix defined by (|4. 15j) . Note that, for i = 1, . . . , n, each 9i, pi 
is a 2 x 2 matrix (since the 9i, pi are eigenfunctions of L, M and L\ M< respectively). 
For integers i,j = 1, . .. ,n, denote by R(i,j) the quasigrammian [11] 



R(i,j) = (-iy 



n(o,p) pt(j) 
e(i) 
















(4.27) 



so that, by once again choosing a trivial vacuum for simplicity, (|4.26p can be expressed as 

S = R(0,0). (4.28) 

As in the quasiwronskian case, we apply the quasideterminant expansion formula (|3.ip . choosing 
the matrices 9i (i = 1, . . . , n) as in (|4.1ip and P = QH^ , where H is a constant square matrix, in 
this case In x 2n, which we assume to be invertible, with denoting the Hermitian conjugate 
of H. Thus 



S 



o(e,p) 







0(6, P) Hipl 



0(6, P) Htf 
ip fo 



0(6, P) 
if) fo 



(4.29) 
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where </>, ^ denote the row vectors (fa, . . . , fan), (fa, • • • > fan) respectively, which gives, by 
comparing the above matrix with (12, 7\i . quasigrammian expressions for q, r, namely 



q = 2a 



n(e,p) 



-2a 



0(9, P) Hfa 



(4.30) 



Thus we have obtained, in (|4,13p . expressions for q, r in terms of quasiwronskians, and in (|4.30p . 
expressions in terms of quasigrammians. We now show how these solutions can be verified by 
direct substitution by firstly explaining the procedure used to determine the derivative of a 
quasideter minant . 



5 Derivatives of a quasideterminant 

We consider a general quasideterminant of the form 

A B 



C 



D 



(5.1) 



where A, B, C and D are matrices of size 2n x 2n, 2n x 2, 2 x 2n and 2x2 respectively. If A is 
a Grammian-like matrix with derivative 



A! = Y i E i F i , 



(5.2) 



i=i 



where Ei (Fi) are column (row) vectors of comparable lengths, it can be shown that (see Ap- 
pendix) 

A R A T3> JL A TP. A U 

(5.3) 



A 


B 




A 


B' 


k 

+£ 

i=l 


A 


Ei 




A 


B 


a 


D' 


+ 


C 





C 







Ft 






where '0' denotes the 2x2 matrix ( q o ) • ^ ^ does not have a Grammian-like structure, we find 
that 



A B 



a d' 



n-l 


A 


fk 


e-k 


+ E 


C 








k=0 









.4 

{A 2 ' 



B 



n—2k\l 



fg2n-2k-ly 
( j^2n—2ky 



(5.4) 



where A k denotes the k th row of A. The formulae (|5.3p and (|5.4p can be utilised to obtain ex- 
pressions for the derivatives of the quasideterminants Q(i,j) (detailed in the Appendix), namely 



Q(i,j)y = Q(i + l,i) - Q(i,3 + 1) + Q(h 0)Q(0, j), (5.5a) 
Q(i,j)x = -a (JQ(i + 1, j) - Q(i, j + 1)J + Q(i,0)JQ(0,i)) , (5.5b) 
g(i, j)t = i (<W + 2,i) - j + 2)J + Q(i, l)JQ(p,j) + Q(t,0) JQ(1, j)) . (5.5c) 

It turns out that the derivatives of R(i,j) match exactly with those of Q(i,j), hence subsequent 
calculations to verify the quasiwronskian solution of the ncDS equations will also be valid in the 
quasigrammian case, meaning that we need only verify one case. 
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6 Direct verification of quasiwronskian and quasigram- 
mian solutions 

We now show that 

5 = Q(0,0) and S = P(0,0) (6.1) 

are solutions of the ncDS system (|2.5al -d). where S is the 2x2 matrix given by (|2.7p and 
A = [J, aS]. Using the derivatives of Q(i,j) obtained in section [5j we have, on setting i = j = 0, 



S y = Q(0, 0) y = Q(l, 0) - Q(0, 1) + Q(0, 0) 2 , 
& = Q(0, 0) x = <t (JQ(1, 0) - Q(0, 1) J + Q(0, 0) JQ(0, 0)) , 
S t = Q(0, 0) f = i ( JQ(2, 0) - Q(0, 2) J + Q(0, 1) JQ(0, 0) + Q(0, 0) JQ(1, 0)) , 
S ra = Q(2, 0) + Q(0, 2) - 2Q(1, 1) - Q(0, 1)Q(0, 0) + Q(0, 0)Q(1, 0) 

+ 2{Q(1, 0)Q(0, 0) - Q(0, 0)Q(0, 1)} + 2Q(0, 0) 3 , 
5 ra = a 2 (Q(2, 0) + Q(0, 2) - 2 JQ(1, 1) J - Q(0, 1)Q(0, 0) + Q(0, 0)Q(1, 0) 
- 2{Q(0, 0) JQ(0, 1) J - JQ(1, 0) JQ(0, 0)} + 2Q(0, 0) JQ(0, 0) JQ(0, 0)) . 

Substituting the above in <\2.9\i . all terms cancel exactly and thus the quasiwronskian solution 
S = Q(0, 0) is verified. As mentioned previously, we obtain the same derivative formulae whether 
we use the quasiwronskian or quasigrammian formulation, and hence the above calculation also 
confirms the validity of the quasigrammian solution S = R(0,0). 



(6.2a) 
(6.2b) 
(6.2c) 

(6.2d) 
(6.2e) 



7 Dromion solutions 



To obtain dromion solutions of the system of ncDS equations, we use the quasigrammian solution 
S = i?(0, 0) rather than the quasiwronskian solution since verification of reality conditions is 
simpler. 



7.1 (n, n)-dromion solution - noncommutative case 



We modify the approach of [10], where dromion solutions of a system of commutative DS equa- 
tions were determined. We consider the ncDS system (|2.5al -d) and, by specifying that certain 
parameters in the quasigrammian are of matrix rather than scalar form, we are able to obtain 
dromion solutions valid in the noncommutative case. Due to the complexity of this solution 
compared to the scalar case considered in jlOj . we look in some detail only at the simplest cases 
of the (1,1)- and (2, 2)-dromion solutions. We do however verify reality for the general case. 
Note here that to obtain dromion solutions, we consider the DSI case, and hence choose a = — 1. 



Recall the expressions for q, r obtained in terms of quasigrammians in (|4.30j) . namely 



o 



12(0, P) Htf 



(7.1) 
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where </>, tp denote the row vectors {<t>i, ■ ■ ■ ,<t>2n), i^ii ■ ■ ■ > V^n) respectively and H = (hij) is 
a constant square invertible matrix, with t denoting conjugate transpose. By once again con- 
sidering the dispersion relations for the system, we are able to choose expressions for (ft, tp 
corresponding to dromion solutions. From (|A.9j) and the definition of = (9±, . . . , 9 n ), where 9i 
is given by (|4.1ip . it follows that <j>, ip satisfy the relations 

{4>j)x = (</>j) y , (4>j)t = K<t>j)yy, ( 7 - 2a ) 

(tpj)x = -(i>j)y, (if>j)t = -'Mj)yy (7.2b) 

(Since the dispersion relations for P are the same as those for 9 when a = -1 (see (TO]) . (f02]) ). 
considering P rather than B and recalling that P = QH^ will give the same relations (|7,2D ). 

So far we have not specified the nature of the noncommutativity we are considering. One 
of the most straightforward cases to consider is to express our fields q and r as 2 x 2 matrices. 
Thus for dromion solutions in the noncommutative case, we choose 

4>j = otjh, (7.3a) 
^ = 0jh, (7.3b) 

where I2 denotes the 2x2 identity matrix, and ay, (3j the exponentials [10^ I13j 

otj = ex.p(pjx + ipjt + pjy + ctj ) , (7.4a) 
/3j = exp(qjx - iqjt - qjy + /3 jo ), (7.4b) 

for j = 1, . . . , 2n, suitable phase constants aj , (3j and constants pj, qj, whose real parts are 
taken to be positive in order to give the correct asymptotic behaviour. The matrix H can be 
assumed to have unit diagonal since we are free to choose the phase constants a JO , /3j Q arbitrarily 
|10j . Using the coordinate transformation X = x + y, Y = —(x — y), we have 

a,j = exp(pjX + ipjt + a j0 ), (7.5a) 
Pj = exp(-^y - iqjt + (3 jo ), (7.5b) 

so that aj = ctj(X,t), j3j = (3j(Y,t). 

We now choose to simplify our notation so that we are working with only <f)\ , . . . , (f> n and 
ipi, . . . , ip n by relabeling (pj as <j> j+i for odd j (i.e. j = 1, 3, . . . , 2n — 1) and setting <j>j = 
for even j ( j = 0, 2, . . . , 2n), and similarly relabeling ipj as for even j and setting tpj = for 
odd j, so that 9j = diag(4>j,tpj) (j = 1, . . . , n) and 

<t> = (0i 02 ... 4>n 0) , (7.6a) 

V> = (0 Vl ^2 ••• lP n ), (7.6b) 

where each <j>j, ipj is a 2 x 2 matrix as defined in (|7.3|) above. Thus, for n = 1, g (which we 
henceforth denote by q 1 for the (1, l)-dromion case and q n for the (n, n)-dromion case) can be 
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expressed in quasigrammian form as 



0(9, p) 



H 



o 



ft* 









0t 



Q 1 



































(7.7) 



where H = (hij) is a constant invertible 4x4 matrix. Applying the quasideterminant expansion 
formula (13, lj) allows us to express g 1 as a 2 x 2 matrix, where each entry is a quasigrammian, 
namely 



0(9, p) 

ai 

0(9, P) 
ai 



h 43 (3t 



0(9, p) 



ai 



0(9, P) 



ai 



/l!4/?l* 











^21 



912 
922 



say. 



(7.8) 



We consider each quasigrammian in turn, but as an example we will look at the quasigrammian 
q\ x . We apply f|3.5 j) to express as a ratio of determinants, namely 

/113# 

0(9, P) : 

^43^1* 

ai „G 
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say. 



(7 9) 

|0(9,P)| F ' " v ' ; 

(We have introduced the notation and <7™ w (u, to = 1, 2) to emphasise that we are considering 
the (u,u>) th entry of the expansion of q n in the (n, n)-dromion case). Although (|3.5|) is valid 
only in the commutative case, our assumption here is that the variables q, r in the system of DS 
equations, and also the parameters 4>j, ipj, are noncommutative. The exponentials ay, (3j given 
by (|7,5p are clearly commutative by definition, hence we are free to use the result (|3.5p . 
By expanding the quasigrammian r in (17. ip in a similar way and extracting the (1, l) th entry, 
we obtain 



r 1 



0(9, P) 



ft 



h\ia\ 




K 
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|0(9,P)| 



F 



say, 



(7.10) 
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with similar results for q\ 2 , r\ 2 etc. 



7.2 Reality conditions 

In the (n, n)-dromion case, we have 



s~<n 

n^VW 

F 

K n 

r) VW 

1 F ' 



(7.11a) 
(7.11b) 



for v,w = 1,2. To verify reality, we must check that r\] m = ±(q™ w )* , i.e. that F is real, and 

(G™ w )* = ±K™ W , where (G™ w )* denotes the complex conjugate of G™ w . 

We do not give details here, however we can show that F can be expressed in the form 



F=\I 4n + H<S>\, 

where H is 4ra x 4n, constant and invertible, and $ is the 4n x An matrix 



(7.12) 



(S 



1^1 





dX 







V 



o 



dX 



I- 



X 

oo 1 







dX 



I- 



X 
oo fn 



*J n dX 




\ 

Jy^ndY 





(7.13) 



remembering that each </>,, ipi (i = 1, . . . , n) is a 2 x 2 matrix and '0' denotes the 2x2 matrix 
(on)- With F in this form, it is straightforward to show that F is real so long as H is a 
Hermitian matrix (see [20]). This condition is also required for (G™ w )* = ±K™ W - in fact, we find 
that r™ w = —(qy W )* so long as H is Hermitian. This agrees with the work of Gilson and Nimmo 
in HDl. 



7.3 (1, l)-dromion solution - matrix case 

We now show computer plots of the (1, l)-dromion solution in the noncommutative case, where 
we choose <j>\ , ipi to be 2 x 2 matrices as in (|7.3p . A suitable choice of the parameters p\ , q\ and 
of the 4x4 Hermitian matrix H allows us to obtain plots of the four quasigrammian solutions 
#11 > Q121 921) #22 as detailed in ()7.8j) . 

We are restricted in our choice of pi, q\ and H in that we require In particular, we derive 

conditions so that F > 0. The determinant F can be expanded in terms of minor matrices of 
H = (hij) = 1,... ,4), giving 

F = 1 + Pi (h 2 H + h\H) e 2 " + Qi (MS + Ml) e" 2 « + PfW 

+ Q 2 M 2 e- 4 « + PxQ! (/l 24 + h\j + + ^13) e 2„-2 € 

+ P 2 Q! (/* 4 + h 3 3 ) e 4 ^ 2 « + PQ 2 (/i 2 + h{) e 2 ^ + PfCfihe 4 *- 4 *, (7.14) 



13 



where Pi 

h rs.. 



l/25K( Pl ), Qi = l/20t(gi), t? = K(pi)(.X - 23( Pl )t), £ = 9\( gi )(Y - 23( qi )t), and 
denotes the minor matrix obtained by removing rows i,j, . . . and columns r,s, . . . of H, 
where i,j,..., r,s,... € {1, 2, 3, 4}. We have used '/i' to indicate that no rows or columns have 
been removed, i.e. h = detH. We can also obtain expressions for each G\, w (v,w = 1,2), for 
instance 

(7.15) 



with ai,/?i defined as in (|7.5p . Similar expansions can be obtained for G\ 2 ,G 2 i and G 22 . Thus, 
it can be seen that for F > 0, we require 9t(pi), £H(gi) and each of the minor matrices in the 
expansion of F to be greater than zero. Suitable choices of p\,q\ and H have been used to obtain 
the dromion plots shown later. In addition to the matrix-valued fields q and r, there are also 
matrix- valued fields A\ and A 2 in the ncDS system (|2.5al d). Plotting the derivatives of these 
fields gives plane waves as follows. 

From (|4.29p . we have an expression for the 2x2 matrix S in terms of quasigrammians. In the 
DSI case A = —\S x — \JS y by (|2.8j) . therefore substituting for S using (|4.29p and equating matrix 
entries gives quasigrammian expressions for A\, A 2 , namely 



A x = -i 
A 2 = -i 



n He 



tp o 



+ i 



4> o 



(7.16a) 
(7.16b) 



We choose <pi, ip± to be 2 X 2 matrices as before so that the boxed expansion element '0' in each 
quasigrammian is the 2x2 matrix ([jjj). Thus, expanding each quasigrammian in the usual 
manner gives a 2 x 2 matrix where each entry is a quasigrammian, and hence we have distinct 
expressions for A\, A 2 corresponding to each of the four dromions q\ t , q\ 2 , q^i-i %2- Considering 
(|2.5cl d). we find that plotting the combination (d x + a d y )A\ for a = — 1 gives a plane wave 
travelling in the X direction, while the combination (d x — crd y )A2 gives a plane wave in the Y 
direction. These, along with the dromions corresponding to each plane wave, have been plotted 
at time t = in figures 1 and 2. 

As can be seen from figure 1, single dromions of differing heights occur in each of the fields 
qili 9i2) 921 and g22- If we were to plot the (1, l)-dromion solution in the commutative (scalar) 
case (that is, if we were to choose q and its complex conjugate to be of scalar rather than 
matrix form), we would obtain only one dromion in the single field q. This dromion and its 
plane waves would have the same basic structure as those above, and thus there would be no 
marked difference in the appearance of the dromions in the commutative and noncommutative 
cases. The main difference between the two situations concerns the number of parameters - a 
far greater number in the noncommutative case gives us more freedom to control the heights of 
the dromions, however some extra care has to be taken in choosing the parameters so that no 
singularities occur in the solution. 



7.4 (2, 2)-dromion solution - matrix case 

In the scalar case [10] , Gilson and Nimmo carried out a detailed asymptotic analysis of their 
(M, A r )-dromion solution, and were able to obtain compact expressions for the phase-shifts and 
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changes in amplitude that occur due to dromion interactions. They then used the results of 
this analysis to study a class of (2, 2)-dromions with scattering- type interaction properties. The 
Hermitian matrix H could be chosen in such a way so that some of the dromions had zero 
amplitude either as t — > — oo or as t — ► +00. 

For the (2, 2)-dromion solution in the matrix case, detailed calculations of this type are more 
complicated due to the large number of terms involved. However, we can adopt the same ap- 
proach to carry out some of the more straightforward calculations. In particular, we obtain plots 
of the situation in which the (1, l) th dromion in each of the solutions q\ x , qf 2 , gfi and Q22 does 
not appear as t — > — 00. These are depicted in figures 3-5. 

To analyse this situation, we focus our attention on g^ and consider G\ x in a frame moving 
with the (1, l) th dromion. We define 

X = X - 23(pi)t, (7.17a) 
Y = Y — 29f(gi)t, (7.17b) 

and consider the limits of Gf 1 as t — > — 00. Let 

m = »(P2) (X - 29f(p2)t) 

= ^(pa) (1 - 2(9(p2) - 3(pi))t) , (7.18a) 

and similarly 

6 = »(<&) (y - 23(<fc)i) 

= 5R(g 2 ) - 2(3(g 2 ) - 3( gi ))t) . (7.18b) 

We choose to order the pi, qi (i = 1, 2) by means of their imaginary parts, so that 9(pi) > 9 (2*2) 
and Q?(gi) < ^fe)- Thus, as t — ► — 00, 772 — » —00 and £2 - * +00. It can easily be shown 
that 772, — £2 determine the real parts of the exponents in 02, $2 respectively, where 02, P2 are 
defined as in (|7.5p . so that, as t — > —00, 02, /?2 — > (and hence a 2 , f3 2 — ► also). Therefore, by 
setting a2 , «2 , > /?2 — ► in and expanding the resulting determinant, we obtain a compact 
expression for G\ x as t — *■ — 00, namely 

r 2 — n R* (h 1245678 p 1 145678 „2t? n . 125678 -2£ , p n , 15678 2/y-2A f 7 i Q \ 

G n - -ai^ ^2345678 ~ - P l /l 345678 e ~ Vl«235678 e + l /i 35678 e ) • I 7 - 19 ) 

Similar expressions can be obtained for the other three determinants G\ 2 , G\\ and G22 by consid- 
ering an extension of (|7.8p to the (2, 2)-dromion case and interchanging columns appropriately: 
for example, we interchange columns 3 and 4, and 7 and 8, in the expansion of G\ x as t — > —00 
to obtain an analogous expression for G\ 2 . Thus we have, as t — ► —00, compact expressions for 



16 



the minors of H governing the (1, l) th dromion in each of Gf 1 , G\ 2 , ^22' namely 



'^2 



+ QlT 



r 2 

U 12 

r 2 

U 22 



h\3 

h^3 



/144 



^23 ^24 
/l43 ^44 




+ PlT 



hi2 

h>22 



hi3 

/l23 



^•13 


h u 




^33 


/134 




^23 


/l24 




^33 


h 3 i 


J 



e'^+PxQxT 



h2i 
h\2 h 13 

h22 ^23 
/l42 /143 



^13 
h23 



h\2 
h22 

/ill 
/l21 



/l!4 
/l24 

/l!4 
^24 



/ill 

/l21 
/l 4 l 



/ll3 
/l23 
/l43 



/l!4 

/l24 
/144 

/l!4 

/l24 
/144 



h\2 


^13 


hu 




h 2 2 


^23 


h 2 4 




h32 


1^33 


^34 




/ill 


hn 


hu 




/i21 


h 2 3 


^24 




^31 


h-33 


^34 


) 



Te 2 V -2Z 



(7.20) 



where T = diag(— 1,1). Since we have written out each minor matrix explicitly, rather than 
using the abbreviated notation as in (|7.19p above, it can easily be seen that, by setting each 
of hi3, hu, /i23j ^24 equal to zero, the (l,l) th dromion in each of <Zii , <7i2 j an< ^ ^22 w * n van " 
ish as t — > —00. This is shown in figures 3-5 below, where we have chosen pi = q2 = \ + i, 
P2 = Ox = 2 ~ an d appropriate values of hij (i, j = 1, . . . , 8). (Note that in the plots of q\ x and 
qf 2 in figure 3, two of the dromions have very small amplitude). We have also shown, in figure 
6, a close-up of one of the dromion interactions at t = 0. 



8 Conclusions 

In this paper we have derived a noncommutative version of the Davey-Stewartson equations 
and verified their quasiwronskian and quasigrammian solutions by direct substitution. The 
quasigrammian solution has then been used to obtain dromion solutions in the matrix case, 
which, if we were to consider these solutions in the scalar case, agree with the results of Gilson 
and Nimmo in pi)]. We have obtained plots of the (1, l)-dromion solution and its plane waves, 
choosing the entries of the Hermitian matrix H in such a way that our solution is always well- 
defined. In the (2, 2)-dromion case, some of the more straightforward asymptotic calculations 
have been carried out, enabling us to obtain plots of the situation with one dromion vanishing 
as t — ► —00. 
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Figure 5: (2, 2)-dromion plots at t = 10. 




Appendix 

Here we prove the results of section [5j Consider a general quasideterminant of the form given 
by (|5,ip . Using the product rule for derivatives, 

E' = D' - C'A~ l B + CA^A'A^B - CA~ l B'. (A.l) 
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We modify slightly the approach of [TTJ and find that, if A is a Grammian-like matrix with 
derivative 



M = Y i E i F i , 



(A.2) 



i=l 



where Ei (Fi) are column (row) vectors of comparable lengths, then 

k 

E' = D' - CA~ l B + ^{CA^E^FiA^B) - CA~ X B' 



A B 

a 



D' 



i=l 

A B' 



C 



i=i 



A Ei 



C 



A B 

F. 



(A.3) 



Note here that '0' denotes the 2x2 matrix ({jo)- If> on the other hand, the matrix A does 
not have a Grammian-like structure, we can once again write the derivative H' as a product of 
quasideterminants as above by inserting the 2n x 2n identity matrix in the form 

n-l 



J = X)(/*e fc )(/i 



kdk) 



(A.4) 



fc=0 



with efe, fk defined as before, so that (fk &k) denotes the 2n x 2 matrix with the (2n — 2k) th and 
(2n — 2k — l) th entries equal to 1 and every other entry 0. Then we find that 



A B 



a d' 



n-l 


A 


fk 


efc 




C 








fc=0 









A 

2n-2fc-ly 



B 



rj^in—zk—i\ 

^j^2n—2ky 



(B 



2n-2fe-ly 
2n-2fe\/ 



(A.5) 



where ^4 fc denotes the k th row of A. It is also possible to obtain a column version of the derivative 
formula by inserting the identity in a different position. We now use the formulae (|A.3P and 
(1A,5|) to derive expressions for the derivatives of the quasideterminants Q(i,j), R(i,j). 



Consider the quasiwronskian Q(i,j) defined in section B~2l namely 

8 

Q(^j) = Q(n+i) 



/i 


e i 















(A.6) 



Calculation of the derivatives of Q(i,j) requires knowledge of the following result [TT], that for 
arbitrarily large n, 



Q{i,j) 



h i+j + 1 = 0, 

i < or j < and i + j + 1^0. 



o o 

. o o 



(A.7) 



We utilise the dispersion relations for the ncDS system (|2.5al -d). found by considering the Lax 
pairs (|2.1al b) in the trivial vacuum case, giving, for 6 an eigenfunction of L, M, 



Ox — —crJOy, 

9 t = iJOyy, 



(A.8a) 
(A.8b) 
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and since = (9\, . . . , O n ) , it follows that 

Q x = -aJQ y , (A.9a) 

& t = Ue yy . (A.9b) 

Thus, using (|A.5|) . we have 

n-l 

Q(i, j) y = Q(i + l, i) + Yl k )Q(~ k > j)> ( A - 10a ) 

A;=0 

/ n-l \ 



Q(i,i)x = + i,j) + J2Q(^ k ) J Qi- k ^) » ( A - 10b ) 

V fc=0 / 

Q(i,j0t = iMQ(t + 2,j0+£Q(* ) *)^(-*,j)) • ( A -!0c) 

V fc=0 / 

These can be simplified using (|A.7|) . leaving the derivatives as given in (|5.5p . 



We can apply a similar procedure to determine the derivatives of the quasigrammian R(i,j) 
defined in section 14.41 The dispersion relations are found by considering the adjoint Lax pairs 
(I4.14al b) in the trivial vacuum case, giving, for p an eigenfunction of L', M', 

p x = ~-Jpy, (A. 11a) 

Pt = lJpyy, (A.llb) 

and since P = (pi, . . . , p n ), we have 

P x = -\j p y, ( A - 12 a) 

P t = UPyy. (A. 12b) 

We also recall that from our construction of the binary Darboux transformation in section 14.3} 
the potential Q(<f>,ip) satisfies the relations (|4.15al c). from which it follows that 

n(e, p) y = pt , (A.i3a) 

fi(9, P) x = -aP^JO, (A.13b) 

n(e, p) t = i (pt J0 (i) _ pW)jq\ , (A.i3c) 

where ^ denotes the k th y-derivative. Thus, using ()A.3p . we calculate the derivatives of R(i,j) 
and find that they are identical to those of Q(i,j) in (|5.5p . 
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